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ABSTRACT 

Imaging of the Sunyaev-ZePdovich (SZ) effect in galaxy clusters combined with cluster plasma x-ray 
diagnostics can measure the cosmic distance scale to high redshift. However, projecting the inverse- 
Compton scattering and x-ray emission along the cluster line-of-sight introduces systematic errors in the 
Hubble constant, Hq, because the true shape of the cluster is not known. In this paper I present a study 
of the systematic errors in the value of Hq, as determined by the x-ray and SZ properties of theoretical 
samples of triaxial isothermal "beta" model clusters, caused by projection effects and observer orientation 
relative to the model clusters' principal axes. I calculate three estimates for Hq for each cluster, based 
on their large and small apparent angular core radii, and their arithmetic mean. I demonstrate that the 
estimates for Hq for a sample of 25 clusters have limited systematic error: the 99.7% confidence intervals 
for the mean estimated Hq analyzing the clusters using either their large or mean angular core radius 
are within ~ 14% of the "true" (assumed) value of Hq (and enclose it) , for a triaxial beta model cluster 
sample possessing a distribution of apparent x-ray cluster ellipticities consistent with that of observed 
x-ray clusters. This limit on the systematic error in Hq caused by cluster shape assumes that each sample 
beta model cluster has fixed shape; deviations from constant shape within the clusters may introduce 
additional uncertainty or bias into this result. 

Subject headings: galaxies: clusters : general, cosmology: distance scale 



1. INTRODUCTION 

There has been a substantial effort to detect the 
Sunyaev-Zel'dovich (SZ) effect from galaxy clusters (Sun- 
yaev & Zel'dovich 1972) and to analyze its distortion of 
the cosmic microwave background radiation (CMB) in con- 
junction with cluster x-ray properties to derive the clus- 
ter cosmological angular-diameter distance and thus esti- 
mates of the cosmological parameters Ho and qo (Gunn 

1978, Silk & White 1978, Cavaliere, Danese, & DeZotti 

1979, Birkinshaw 1979; see also Birkinshaw 1998, and ref- 
erences therein). This method provides a distance deter- 
mination for the cluster that is independent of the "cosmic 
distance ladder" of Cepheid variables or supernova?, and 
is potentially effective for clusters at high redshift (z ~ 1). 
Centimeter- wavelength interfcrometry optimized for imag- 
ing the SZ effect from galaxy clusters has been recently de- 
veloped (Carlstrom, Joy, and Grego 1996; Grainge et al. 
1996;). This allows high-resolution x-ray and radio images 
of clusters to be analyzed simultaneously. The results of 
fits of both the x-ray and radio images to simple cluster- 
plasma models will yield improved estimates of Hq, and 
systematic errors in the measured value of Hq are likely to 
be a significant limit to its accuracy. 

Sources of systematic errors in the SZ-detcrmined Hq 
and qq can originate from the assumptions made in mod- 
eling the cluster plasma: ignorance of the cluster plasma's 
true three-dimensional distribution and inadequate treat- 
ment of the physical state of the cluster plasma. Radio 
and x-ray images only provide the projected x-ray surface 
brightness and CMB decrement. For the analysis to pro- 
ceed some assumption must be made about the cluster 



size along the line of sight; e.g. , one assumes that clus- 
ter has spherical or ellipsoidal symmetry. The modeling 
of physical state of cluster plasma for SZ analysis has gen- 
erally assumed that the plasma was of a single phase and 
temperature, using the somewhat ad hoc "beta" model for 
electron density, n e (r) oc (l+r 2 /r 2 )~ 3 ' 3 / 2 , where r c is clus- 
ter's "core radius", within which the density is relatively 
flat. The beta model can be argued as a possible distri- 
bution for the plasma in a dynamically relaxed isother- 
mal cluster in hydrostatic equilibrium (e.g. , Cavaliere & 
Frusco-Femiano 1978; Sarazin 1986), but its usefulness is 
more empirical; many x-ray images of clusters fit a beta 
model reasonably well (Mohr et al. 1999; see §4). Also, 
studies of the distribution of SZ systematic errors caused 
by cluster shape and orientation (and other effects) based 
on the results of a large ensemble of numerically simu- 
lated clusters have yet to be completed; current results 
are for a small set of simulated clusters (see §4). Thus, 
three-dimensional "toy model" estimates for the effects of 
cluster shape are a useful first step in estimating these er- 
rors, and can help identify the physical sources of bias and 
scatter in Hq estimates from simulated clusters. 

In this paper I study the systematic errors in the value 
of Hq, measured by SZ and x-ray observations, caused 
by effects of cluster shape. This study consists of two 
parts. First, I create theoretical galaxy cluster samples, 
where each cluster's plasma distribution follows a triaxial 
isothermal beta-model (§2), possessing three independent 
core radii. I use the triaxial beta model because it is a sim- 
ple three-dimensional generalization of the spherical or el- 
lipsoidal beta models (commonly used in SZ analysis) that 
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demonstrates the effects of shape and orientation on the 
uncertainties in Hq determined by SZ observations. The 
triaxial beta model also produces simple analytical func- 
tions for the CMB decrement and x-ray surface brightness 
so results for large samples of clusters can be easily calcu- 
lated. 

I create numerical distributions clusters by uniformly 
and independently sampling the plasma core radii, con- 
straining them by a minimum ratio between any two core 
radii of a sample cluster. These samples are uniform in 
the plane of allowed cluster oblateness and ellipticity. The 
clusters are placed in the sky with a random orientation 
to our line of sight. I identify a cluster sample with an 
optimum asphericity that has a distribution of apparent 
cluster cllipticities that is consistent with that of observed 
x-ray clusters (Mohr et at 1995; see §3). 

Second, I analyze the clusters of the theoretical sample 
to determine their distance as if they were either spheri- 
cal or an ellipsoid of rotation, as in done in observational 
analysis {e.g. Hughes & Birkinshaw 1998). An important 
unknown quantity is an ellipsoidal cluster's inclination an- 
gle i; the estimated value for Hq will vary greatly with i. 
However, since our theoretical clusters are actually three- 
dimensional, specifying a single inclination angle is artif- 
ical. Therefore, I analyze each cluster very simply as if 
its inclination angle i — 90°, i.e. , that the core radii for 
the clusters are not altered by projection effects, and then 
study the distribution of the estimates for Ho for a large 
number of sample clusters. The apparent shape of a sam- 
ple cluster's x-ray surface brightness will be elliptical, with 
a large and small angular core radius, 9 + > 8~ . I calculate 
two different estimates of Hq which are proportional to ei- 
ther + or 9~ , designated Hq and Hq . I also calculate an 
estimate Hq VS , by using the arithmetic average of 9 + and 

I find that the sample means of the estimates Hq and 

Hq VS fall within ~ 5% of Hq for the optimal sample. The 

sample distribution of Hq shows greatest bias, with a 

mean for Hq that underestimates Hq by ~ 14% for my 
optimal cluster sample. 

As a predictor for SZ observations, I also calculate es- 
timates for Ho averaged for 1000 realizations of a sample 
of 25 clusters. I find that the systematic errors caused by 
cluster shape are limited: the 99.7% confidence intervals 
for Hq and Hq VS include the assumed value of Hq for my 
optimal cluster sample, and do not extend beyond ~ 14% 
from Hq. The 99.7% confidence interval for Hq does not 
include Hq, indicating that it may not be a useful param- 
eter for distance estimation. 

The structure of this paper is as follows. In section §2 
I describe the triaxial beta model for the cluster plasma, 
and describe the analytic expressions for their CMB decre- 
ment and x-ray surface brightness. I also describe the con- 
struction of samples of theoretical clusters, distinguished 
by their degree of triaxiality, and describe the manner in 
which I analyze the apparent clusters to determine values 
for Hq. I present our results in §3, followed by a summary 
and discussion - noting some of the limitations of this beta 
model based analysis - in §4. 

2. METHOD 



2.1. Triaxial beta model clusters 

The distribution of cluster plasma is described by an 
isothermal "beta" model. The electron density at a posi- 
tion within the cluster r = xixi + £2X2 + £3X3, measured 
in the observer's coordinates, is given by 



; (r) = n e0 (l + r-M-r) 



(1) 



where the matrix M describes a cluster's shape and ori- 
entation of its principal axes to the observer, and (3 is an 
exponent with the nominal range | < /3 < 1. The maps of 
x-ray surface brightness Sx and cosmic-microwave back- 
ground decrement 5T r in sky angular coordinates tp) 
(measured from the cluster center r = 0) are given by in- 
tegrals of the x-ray emissivity and electron pressure over 
the line-of-sight (defined here as x%) though the cluster; 
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T e (r)n e (v)dl. (3) 



Here T e is the electron temperature, hereafter assumed 
to be constant, Ax(T e ) is the plasma emission function 
over a prescribed x-ray bandwidth at temperature T e , and 
z is the cluster redshift. 

For a triaxial isothermal beta model plasma described 
by equation ([!]) , then integrating equation (^|) by choosing 
dl = dx\ gives Sx to be 



Sx{#,<p) = 



B(§,3/?-§) nl e A x (T e ) 



J eff 



4?r 



(i + ^) 3 x($^Y p - 



(4) 



where x($i f) is a quadratic function of the sky angular 
coordinates d and ip, describing elliptical isophotcs. Along 
the line of sight of the center x(0, 0) = 1- The quantity 
B(q, r) is the beta function. The quantity L e ff is an effec- 
tive column length for the plasma along the line of sight 
through the cluster: 
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The quantities (r c i, r C 2, 7*03) are the cluster core radii, 
and (ai, 0.2,013) are the rotation angles of the cluster prin- 
cipal axes relative to the observer. The coefficients in the 
quadratic function x($i f>) are also functions of the cluster 
core radii and its orientation. 

By integrating equation (^) in exactly the same manner, 
ST r can be shown to be 
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The value of L e ff determined from observations by re- 
lating Sx and ST r observed for the cluster. For example, 
L e f f can be determined by the values of ST r and Sx mea- 
sured at the cluster's center: 

£?(! ,30-1) 1 Ax{Te) / meC 2y 5T 2^ Q) 

eff (1 + z) 3 47ra| \k B Tj S x (0,0)- 

(7) 

The cluster's cosmological angular diameter distance 
Dg(z; Hq, go) is then inferred by equating the measured 
L e ff to that derived from a model for the cluster. Recent 
efforts to fit the cluster Sx (Hughes and Birkinshaw 1998) 
assume that the cluster is either an oblate or prolate ellip- 
soid in shape as well as isothermal; this assumption about 
shape is reasonable when no information can be known 
about the structure of the cluster along the line of sight. 
However, the dependence of £ e // on the cluster's appar- 
ent major and minor axes for an ellipsoidal beta model is 
different than that for a triaxial cluster. 

f 2 ■ • 2 ■ ^ -1/2 

ellipsoidal I COS % Sin I 

L eff = D e{z;H ,q )^-^— + -^-j , (8) 

where 9 C \ and 9 c2 are an ellipsoidal clusters' angular 
axes; 6> cl > 9 c2 describes an oblate ellipsoid and 9 cl < 9 c2 
describes a prolate ellipsoid. The quantity i is the inclina- 
tion angle of the symmetry axis. 

If the cluster is triaxial, the observed L e f f will be equal 
to that of equation (^|). However, analyzing the cluster 
as an apparent ellipsoid with equation (@) will produce a 
systematic error in the inferred value of Dg, and thus in 
cosmological parameters Hq and qo. 

2.2. The theoretical cluster samples 

I generate triaxial beta-model clusters choosing a set of 
core radii (r c i, r c2 , r 3) from a random uniform distribu- 
tion for the ratio of two of the core radii, r c2 and r c3 , 
with respect to r c \ . Both r c2 and r C 3 are assumed to be a 
random fraction of the length of r c i, but bounded below 
by a minimum value. This minimum value is not known 
a priori but can be chosen to optimize the observed el- 
lipticity of x-ray clusters (see below). Clusters are only 
distinguished by their core radii; I do not create a dis- 
tribution for the clusters' (3- values nor any other quantity 
except core radii. Spherical beta model fits to real clusters 
appear to exhibit correlation between core radius and the 
value of (3 (Neumann and Arnaud 1999). However, the ob- 
served correlation for beta model clusters is not convolved 
by cluster shape projection. This is shown in the expres- 
sions for the x-ray surface brightness Sx-, equation (|J), 
and SZ CMB decrement, ST r , equation (g). The profile 
exponents for both of these quantities are not functions of 
the individual core radii nor of rotation angles; nor are the 
major and minor axes of the observed elliptical cluster (9 + 
and 0~ , see §2.3) functions of /3; they are only functions 
of the sample cluster's core radii and the rotation angles. 
Also, I am not considering a distribution of the magnitude 
of the cluster core radii, but the distribution of the cluster 
triaxiality (see §2.3). 

I rule out bias in the cluster samples toward a net oblate- 
ness or prolateness by checking for uniform sampling in the 



ellipticity-prolateness (E, P) plane, given by Thomas et al. 
(1998) as 



P _ 1 r c2( r cl r c3) /q.\ 

~ 2 r c Vc 2 3 + r 2 cl r 2 c3 + r 2 cl r 2 c2 ' 

and 

p_ 1 'c2'c3 Al cV c3, ' ' cl' c2 qq\ 

~ 2 r 2 2 r 2 c3 + r 2 cl r 2 c3 + r 2 cl r 2 c2 ' 

Strictly prolate and oblate clusters fall onto the lines 
P = —E and P = E respectively, with length determined 
by the lower bound of the ratio between core radii. My 
optimal sample, described below, uniformly covers the al- 
lowed region in the {E, P) plane, which is a triangle pro- 
scribed by the prolate and oblate lines and the line con- 
necting their endpoints. 

How well does the triaxial beta model reproduce the 
observed shapes of x-ray clusters that could be used for 
SZ analysis? A study of 65 Einstein x-ray clusters by 
Mohr et al. (1995) found an emission-weighted mean el- 
lipticity of 0.20 ± 0.12, while McMillan, Kowalski, and 
Ulmer (1989) found a mean ellipticity of 0.24 ±0.14 for 49 
Einstein Abell clusters. In both of these studies clusters 
were included with substantial flattening caused by recent 
merging, or with cooling inflows which can make the clus- 
ter appear more spherical. A more appropriate sample for 
comparison would be one which excludes these effects. If 
I eliminate clusters that are apparent mergers from the 
Mohr et al. sample (8 out of 12 clusters with ellipticities 
of 0.3 or greater with apparent subclustering) and clus- 
ters in which cooling inflows may exist (as measured by 
central cooling times of 10 Gyr or less; an additional 17 
clusters), then the mean ellipticity of the remaining sub- 
set is 0.18 ± 0.09. I find that a triaxial beta model cluster 
sample where the minimum ratio between any two core 
radii to be ~ 0.65 produces a distribution of apparent el- 
lipticities that is consistent with this subset of the Mohr 
et al. sample (figure 1). A Kolmogorov-Smirnov test be- 
tween these samples indicates the ellipticity distributions 
are statistically indistinguishable, with a maximum differ- 
ence between the two cumulative distributions of d = 0.12, 
and a probability that the two samples are drawn from the 
same distribution of ~ 69%. 

2.3. The analysis 

The systematic error in analyzing the clusters arises 
from assuming that an apparent cluster is either a prolate 
or an oblate ellipsoid, when it is in fact triaxial. The ap- 
parent elliptical image of the cluster will have a major and 
minor axis, + and 9~; these are the "angular core radii" 
for the x-ray and SZ images. Using equation (^) and the 
function x($>¥'); ®~ i and L e ff arc determined for a 
given triaxial cluster. The observational analysis proceeds 
as if the observed 9 + and 9~ are that of an ellipsoidal clus- 
ter, inclined to the line of sight by an unknown angle i. 
The apparent cluster distance Dg for an ellipsoidal clus- 
ter is related to 9 + , 9~ , L e ff, and i by deprojecting the 
cluster axes and using equation (|sj) : 
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(prolate) ; 
(oblate). 

Equating L e fj of equation (O) with that for the tri- 
axial cluster, then in general Dg will not be equal to the 
actual distance Dg. This leads to erroneous values for the 
apparent cosmological parameters Hq and go- 

Since the sample clusters are intrinsically triaxial, us- 
ing estimators (|llj) for Hq based on an ellipsoidal cluster 
model with a single inclination angle i is artificial; there 
is no single angle that characterizes the orientation of a 
cluster to the line of sight, unless a sample cluster's core 
radii had been accidentally chosen to be roughly prolate 
or oblate. Therefore, I collapse the dependency on i and 
assume i = 90°, and use only a cluster's observed 9 + and 
6~ . I compose three estimates for Hq for each cluster: 
H+ oc 6+ , Hq oc 6~, and H™ 6 cx \(6 + + 9~). These 
estimates are ordered Hq > Hq VS > Hq . The esti- 
mate H^ is equivalent to assuming the observed cluster 
an oblate ellipsoid, while the estimate Hq is for the clus- 
ter as a prolate ellipsoid, both viewed as if the axis of 
rotation were in the plane of the sky. I study the distri- 
bution of these estimators for Hq for a sample of triax- 
ial clusters drawn from the distribution described in §2.2. 
The cosmological parameter go is fixed to be zero, so that 
Dg = cHq 1 z(1 + \z)/{\ + z) 2 . As mentioned in §2.2, I 
am sampling clusters only by a distribtution in their tri- 
axiality, and I do not use the magnitude of the core radii. 
Thus the estimators Hq are determined with respect to an 
assumed value of Hq . 

3. RESULTS 

Figure 1 shows the distributions of the values of Hq 
and Hq inferred from the cluster apparent angular axes, a 
scatterplot of these estimates (one point per cluster), and 
the distribution of apparent cluster ellipticities, for our 
optimal theoretical beta model sample (with a minimum 
ratio of core radii of 0.65). Figure 2 shows the distribution 
of of Hq VS for this cluster sample. 

The distributions of Hq has sample mean that is within 
~ 8% of the assumed value of Hq. The distribution of Hq 
is more biased {i.e. , lower) in sample mean value, which 
is expected since for every sample cluster the value of Hq 
is constrained to be lower in value than Hq . 

I also studied the distribution of the Hq estimates for 
clusters samples with greater asphericity. For a cluster 
sample with a minimum ratio of core radii of 0.5, the es- 
timate distributions broaden significantly, and the sample 
deviations nearly double. The distribution of Hq VS also 
exhibit broadening with greater asphericity, but still has 
a mean value within ~ 5% of Hq. All of the means are 
relatively insensitive to the samples' degree of asphericity. 

What then is the expected systematic error in the mea- 
sured Hq caused by cluster shape for a practical-sized sam- 
ple of clusters? The 99.7% confidence intervals for the 
mean values Hq, Hq , and Hq V& , based on 1000 realiza- 
tions of a 25-cluster sample are summarized in Table 1. 



The confidence intervals for Hq and Hq VS include the as- 
sumed value of Hq for the optimal cluster sample, and even 
for samples with even greater asphericity. For the optimal 
cluster sample the confidence intervals for Hq and Hq VS 
do not extend beyond ~ 14% from Hq. The 99.7% confi- 
dence interval for Hq does not include the assumed value 
offlo- 



4. SUMMARY AND DISCUSSION 

The high-resolution imaging of the SZ effect in galaxy 
clusters in combination with cluster plasma x-ray diagnos- 
tics is a powerful technique for measuring the cosmic dis- 
tance scale. This method is sensitive to the projection of 
the cluster's inverse-Compton scattering and x-ray emis- 
sion, which depend on plasma density and temperature 
along the cluster line-of-sight. 

In this paper I have estimated the systematic errors in 
the SZ-determincd Hubble constant caused only by the 
projection effects of cluster shape. I use a triaxial beta 
model to represent the clusters' gas because it is the sim- 
plest generalization of the ubiquitous spherical and ellip- 
soidal beta models that can demonstrate the effects of clus- 
ter shape and orientation on measurements of Hq. The 
triaxial beta model has analytic expressions for the clus- 
ters' CMB decrement and x-ray surface brightness, so that 
the statistics for measured Hq for a very large number of 
clusters are easily computed. 

Ideal clusters for SZ analysis would possess no obvious 
substructure nor evidence of merging, and contain plasma 
at a single temperature without cooling inflows. As most 
observed clusters are not such simple systems, I discuss 
the relevance of my beta model sample predictions for 
errors in Hq caused by cluster shape for a real SZ clus- 
ter survey. First, will the presence of cooling gas alter 
a cluster's SZ properties substantially from a beta model 
description ? A recent analysis of ROSAT x-ray clusters 
indicate that a majority of them contain cooling gas (Peres 
et al. 1998). The clusters' x-ray emission is sensitive to 
the environment of their centers, where gas densities are 
highest, while their SZ effect is relatively more sensitive 
to the lower-density outer regions of the clusters. Cooling 
inflows are confined to the core of the clusters and are be- 
lieved to be subsonic and isobaric, so that the relevant SZ 
cluster property, the column integral of the cluster pres- 
sure, will be largely unaffected with the presence of cooling 
gas in the core. "Reprojection" estimates for Einstein x- 
ray clusters suggest that the SZ effect is enhanced in the 
largest of cooling clusters, with cooling rates of hundreds 
of solar masses per year (White et al. 1997). However, 
most of these clusters have high pressures in their outer 
regions, so the SZ enhancement is likely not to be caused 
by the alteration of the SZ effect within the cooling core, 
but by the presence of higher gas pressure over the ex- 
tended outer (beta model) region. Strong cooling inflow 
clusters have sharply peaked x-ray profiles so that their x- 
ray determined core radii and central plasma densities can 
be skewed. As mentioned in §2.2, in comparing my beta 
model sample to the observed cluster ellipticity distribu- 
tion, I eliminate the stronger cooling inflow clusters in at- 
tempting to avoid this bias in observed shape. For many of 
these systems these issues are academic, as they are radio 
loud and obscure a CMB decrement (Burns 1990). How- 
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Minimum ratio of core radii Hq / Hq H q / Hq 



rravg 



/Ho 



0.65 



0.94 - 1.14 0.78 - 0.97 0.87 - 1.05 



Table 1 

99.7% CONFIDENCE INTERVALS FOR 25-CLUSTER SAMPLE MEAN OF Hq. THE MININUM RATIO BETWEEN CORE RADII IS 0.65 



ever, some of the clusters that have a detected SZ effect 
contain large cooling inflows (Hughes 1997), and observers 
have used the beta model properties of the outer portion 
of the clusters to extract the relevant x-ray properties for 
analysis in conjunction with the SZ effect (e.g. , Myers 
et al. 1997). It is also interesting to note that for many 
clusters - approximately half of ROSAT clusters of an x- 
ray flux-limited sample previously selected by Edge et al. 
(1990), including some that contain cooling inflows - the 
simple beta model can produce a good fit to their x-ray 
profiles (Mohr et al. 1999). 

An important caveat to the error estimates in Hq pro- 
vided by these beta models is that they cannot account 
for the cluster gas distibutution changing shape from the 
core to the outer region. Observations of a few high sig- 
nal to noise ROSAT clusters show ellipticity gradients, 
exhibiting a rougly linear decline in x-ray ellipticity from 
e ~ 0.3 — 0.1 from the clusters' center, over a distance of 
several core radii (Buote and Canizares 1996). This be- 
havior may have a substantial effect on the SZ properties 
of a cluster, altering both the shapes and magnitudes of 
the apparent x-ray and (to an even larger degree) the SZ 
images. While results of cluster ellipticities from a larger 
sample are required to adequately constrain this effect for 
a statistical study, here I illustrate a possible bias in Hq 
that could arise from changing cluster shape using a sim- 
ple model of a cluster with varying ellipticity I consider 
a beta model with the core radii in the elements of M 
in equation (|l|) as functions of coordinates. An ad hoc 
example is 



r 2 
1 cl 



r c i + r 



(i? 2 -r cl ), 



(12) 



used in the beta model of equation ([!]), with where r is the 
distance of the coordinate point from the cluster center. 
This describes a triaxial cluster (using similar expressions 
for r C 2 and r c ^) with core radii of (r c i, r C 2, ?"c3) within the 
cluster core, becoming spherical with core radius R > T c \ 
outside the cluster's core. I refer to this as a cluster with a 
"modified" core radius. Choosing a ~ 4 and r c i/R ~ 0.6 
for a cluster with one core radius modified by equation 
dl2] ) and f3 = | in equation (|l|), produces a decreasing x- 
ray ellipticity from e(R) ~ 0.3 to e(5R) ~ 0.1, when the 
modified core radius lies in the plane of the sky; there is 
similar behavior in ellipticity for two modified cluster core 
radii with one along the line of sight. The presence of gas 
with a more spherical distribution moderates the effect of 
cluster orientation on L e ff determined from equation (Q), 
with L e ff assuming intermediate values within the range 
given by the outer core radius R and the inner core radii 
{r c i,r C 2,r c ^). For example, a cluster with one modified 
core radius, taken along the line of sight and using the 
values for r c \/R ~ 0.6 and a ~ 4, produces an estimate 



Hq that is biased (high) by ~ 20%. This by itself is a 
substantial effect on what otherwise appears as a spheri- 
cal cluster, however in conjunction with orientation, it is 
a significantly lower bias that would have been produced 
by the unmodified oblate cluster observed along its minor 
axis, ~ 67%. I have calculated the x-ray and CMB decre- 
ment images for a small set of prolate or oblate clusters, 
with one or two modified core radii (using the parameters 
from above), observed along the axes and along the line 
x — y — z. For these clusters the estimates Hq and Hq 
are lower (by « 10%), and their difference Hq — Hq sub- 
stantially smaller than their counterparts for clusters with 
unmodified core radii. It is uncertain whether any signifi- 
cant bias to estimates for Hq would be introduced in statis- 
tically analyzing a large set of such type of clusters. These 
estimates are based on using the central values for x-ray 
brightness, CMB decrement, and determining the cluster 
angular size using the apparent x-ray core radii defined 
simply by the ellipse of brightness that is lower than the 
peak by a factor 2 3/3 ~ 1 / 2 with /3 = |. More quantitative 
results will depend on the model details of the distribution 
of gas in transition from core to outer region of the clus- 
ter, and the manner by which models are fit to the data to 
determine parameters (e.g. simultaneous x-ray and CMB 
image fitting), well beyond the scope of this paper. Quali- 
tatively, however, the presence of a changing cluster shape 
can alter the estimates for Hq by softening the effects of 
orientation. A notable effect for this type of cluster is that 
the beta model congruence of the x-ray and CMB images 
(equations (||) and (|6j)) is broken, so that comparison of 
the maps may determine the importance of shape changes 
outside of the cluster core. 

Clusters with recent merging activity cannot be ade- 
quately represented by simple beta models. However, the 
use of such clusters in a SZ survey is likely to be fraut 
with complications. In principle, numerical simulations of 
cluster formation would yield a more "realistic" sample of 
clusters for SZ analysis than my beta model sample, ac- 
counting for the effects of cooling and merging as well as 
for shape projection. However, simulations of cluster for- 
mation can not yet physically reproduce the observed large 
gas cores that are observed in clusters (Metzler and Evrard 
1997; Anninos and Norman 1996). Inagaki et al. (1995) 
used two simulated clusters to study SZ measurement sys- 
tematics caused by plasma temperature gradients, plasma 
dumpiness, cluster peculiar velocity, the finite extent of 
cluster plasma, and cluster shape. They determined that 
effects of asphericity would be limited to an uncertainty of 
~ 10% in Hq by observing several clusters, as I have also 
found. They did not conduct a survey of possible cluster 
shapes; the statistics for estimates of Hq were generated 
by the viewing of the two clusters at many orientations. 
I have constrained the limits of my beta model shapes by 
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checking for its consistency with the observed apparent 
cllipicities of x-ray clusters. Roettiger et al. (1997) fo- 
cused on the systematic errors in an SZ-determined Hq 
observed in seven simulated cluster mergers with strong 
temperature gradients and asphericity. They found that 
these effects could lead to H underestimated by as much 
as 35%, and concluded that two approaches should be used 
in SZ analysis: perform detailed simulations of individual 
clusters where it was indicated that merging was strongly 
affecting the SZ properties, otherwise employ a statistical 
sample of clusters that show no evidence of recent merging 
or dynamical evolution. In this paper I have addressed the 
systematic errors caused by cluster shape and orientation 
that would be present in using this latter approach with 
a modeled optimal SZ cluster sample. What are needed 
now are the statistical results for SZ Hq estimates from a 
large sample of numerically simulated clusters. 

I have created numerical samples of triaxial beta model 
clusters by specifying the minimum ratio between any two 
core radii. I have identified an optimal such sample, with 
the ratio of ~ 0.65, that has a distribution of apparent 
cluster x-ray ellipticities that is consistent with that mea- 
sured from observations of x-ray clusters. I have analyzed 
the cluster samples for their SZ decrement and x-ray sur- 
face brightness, assuming no effects of inclination angle. 
The apparent cluster's large and small angular core ra- 
dius, 8 + and 6~ , yield three estimates of Ho that are pro- 
portional to 9 + , 9~ and \{9 + + 6~). These estimates are 
equivalent to assuming that the cluster is either oblate 
(oc 9 + ) or prolate (oc 9~), with its symmetry axis in the 
plane of the sky, or spherical oc \{9 + + 6~). I have found 



that the estimates Hq and Hq V£ , have means that fall 
within ~ 5% of the assumed value of H for the optimal 
theoretical cluster sample, while the mean of the estimate 
Hq underestimates Hq by ~ 14%. The size of these errors 
caused by cluster shape is similar to that found in a more 
approximate fashion by Hughes and Birkinshaw (1998), 
and discussed recently by Cooray (1998). Other estimates 
of Hq can be devised, for example, a (weighted) geometric 
mean Hq = (iJ+) a (i? -) 1 - Q ; < a < 1 (Van Speybroeck 
and Vikhlinin 1997). These may produce better estimates 
for Hq than the three simple estimates that I have stud- 
ied, but the best choice of Ho estimator may depend on 
the intrinsic shape distribution of clusters. 

I have also determined the confidence intervals for the 
estimates of Ho that would be derived from the SZ and 
x-ray analysis of a 25-cluster sample. Our optimal the- 
oretical cluster sample has 99.7% confidence intervals for 
Hq and Hq V& that are within ~ 14% of Ho, and enclose 
Hq- The confidence intervals for the estimate Hq show 
more deviation and do not enclose Ho, indicating that it 
may not a useful estimator. 
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Fig. 1. — The distributions of Hg inferred from the cluster apparent core radii, compared to the assumed value of Hq, a scatterplot of 
if+ vs. H~ , and the distribution of apparent cluster ellipticities, for the optimal sample of 1000 triaxial beta model clusters This sample 
produces a distribution of apparent cluster ellipticities that is consistent with the subset of the Mohr et al. sample (see text). Distribution 
mean values are indicated by dashed lines (value of sample mean and deviation above line); the dashed line on the scatterplot indicates where 
cluster appears circular. The dotted-linc histogram is that of the subset of the Mohr et al. x-ray clusters described in the text (full scale for 
its frequency is — 20). 
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The distributions of H, 







inferred from the clusters by averaging major and minor axes, compared to the assumed value of Hq, 



for the optimal sample. The distribution mean values are indicated by dashed lines (value of sample means and deviations above lines). 



